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Abstract: We study the medium-induced energy loss —AE (L P ) suffered by a 
heavy quark produced at initial time in a quark-gluon plasma, and escaping the 
plasma after travelling the distance L p . The heavy quark is treated classically, and 
within the same framework —AE (L p ) consistently includes: the loss from stan- 
dard collisional processes, initial bremsstrahlung due to the sudden acceleration 
of the quark, and transition radiation. The radiative loss induced by rescatter- 
ings —AE ra d(L p ) is not included in our study For a ultrarelativistic heavy quark 
with momentum p > 10 GeV, and for a finite plasma with L p < 5 fm, the loss 
— AE (L P ) is strongly suppressed compared to the stationary collisional contribution 
—AE co u(Lp) oc L p . Our results support that —AE ra d is the dominant contribution 
to the heavy quark energy loss (at least for L p < 5 fm), as indeed assumed in most 
of jet-quenching analyses. However they might raise some question concerning the 
RHIC data on large p± electron spectra. 
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1. Introduction 

The attenuation of jets - jet quenching - was initially suggested by Bjorken [1] as a 
possible signature of the quark-gluon plasma (QGP). Since, the suppression of hadron 
spectra at large transverse momentum p±, in ultrarelativistic heavy ion compared 
to proton-proton collisions, has been observed [2, 3] at the Relativistic Heavy Ion 
Collider (RHIC). Until recently this suppression (and other features of jet quenching) 
seemed consistent [4] with the energy loss —AE ra d of the parent parton due to gluon 
radiation induced by its rescatterings [5,6] in the hot or dense medium created in the 
heavy ion collision. 

However, recent experimental data on the p^-spectra of electrons arising from 
heavy flavour decays [7,8], when compared to theory [9], suggest that heavy quarks 
lose more energy than —AE ra d only 1 . It was suggested in Ref. [11] that this 'single 
electron puzzle' might be solved by the collisional energy loss, previously neglected 
in the total parton energy loss. However, one should remember that other sources 
of medium-induced energy loss are usually neglected, namely the Ter-Mikayelian 
(TM) effect and transition radiation. Thus it would be incorrect to invoke collisional 
loss without studying the effect of the latter radiative contributions 2 . Transition 

1 It is stressed in Ref. [10] that this statement might be somewhat premature, since the theoretical 
calculation of heavy quark production suffers from large uncertainties already in proton-proton 
collisions, and also because the contributions to the electron spectra from charm and beauty are 
not yet separated experimentally. 

2 We stress again that the TM effect and transition radiation are radiative contributions distinct 
from that induced by rescatterings —AE rac i- 
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radiation and the TM effect have already been studied, for instance in Refs. [12, 
13] and [14] respectively. In this paper we present a theoretical model allowing to 
treat consistently the zeroth order of the (heavy quark) energy loss in an opacity 
expansion [6], denoted as —AE , which includes the TM effect, transition radiation, 
and collisional energy loss. Our main result is shown in Fig. |^. For L p < 5 fm we find 
that —AE is suppressed compared to the collisional loss — AE co u, indicating that 
—AE rac i should be dominant for such values of the plasma size (contrary to what is 
suggested in [11]), and leaving the 'single electron puzzle' unsolved. 

In section we present our theoretical model, which is inspired from that of 
Ref. [15], in the case of an infinite medium. We assume the heavy quark to be 
produced in a hard subprocess at initial time, in a static QGP of high temper- 
ature T and small coupling g <C 1. The latter hypothesis implies the hierarchy 
1/T <C ru <C A, where 1/T is the average distance between two constituents of the 
QGP, r D = l/m D ~ l/(gT) is the Debye radius and A ~ l/(g 2 T) is the mean free 
path of the heavy quark. Under these hypotheses, we can describe the QGP via its 
collective response to the current [16], where the (longitudinal and transverse) dielec- 
tric functions are obtained from the gluon polarization tensor. In [15] we studied the 
mechanical work W undergone by the heavy quark when travelling in the (infinite) 
plasma, and used the hard thermal loop (HTL) [17] approximation for the dielectric 
functions. In section ^ we show that to a very good accuracy, the model of [15] can be 



simplified by assuming an isotropic plasma described by the dielectric function (|2] 



In fact, in the ultrarelativistic limit j — E/M— l/yl — v 2 ^ oo we are considering 
(M, E and v being the heavy quark mass, energy and velocity), the model ( |2.S| ) 
reproduces exactly the results of [15] for the work W. At finite but large 7 > 1, the 
differences between the HTL approximation and the simple model ( |2.8| ) are numer- 
ically negligible, as shown in Fig. [I]. We end section |2| by a critical discussion. We 
indeed stress that in the academic case of an infinite QGP, the mechanical work W 
calculated in [15] should not be interpreted as an observable energy loss (contrary to 
what was done in [15]). This is due to the fact that the quark asymptotic (t = +00) 
states are different in vacuum and in an infinite QGP [18]. 

In section [| we consider the more realistic case where the quark is produced in 
a hard subprocess in a finite size (but still static) QGP, and travels the distance 
L p before escaping the medium. In this case the "asymptotic" quark states are the 
same in vacuum and in the presence of the medium, and it is legitimate to interpret 
the work — W as the (medium-induced) quark energy loss, provided one includes 
transition radiation. We thus study the heavy quark energy loss — W = —AE(L, L p ), 
where L is the total distance travelled by the (deconfined) quark, related to the 
quark hadronization time as L = wthad ~ ^7/Aqcd- We assume the heavy quark 
to be relativistic enough so that it hadronizes outside the medium, thad — L ^> L p . 
We also consider a high temperature QGP, with a Debye mass ~ gT ^> Aqcd- 
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The hadronization time being larger 3 than the other relevant length scales of the 
problem L p and 1/mo, the loss —AE(L = oo,L p ) = —AEq(L p ) will be the main 
quantity of interest, affecting the final hadron (or decay electron) pj_-spectrum. The 
calculation of transition radiation fields in section |3| is done using the simple model 
( |2.8| ) motivated in section 0, and our final result for the zeroth order energy loss 
— AE (L P ) is shown in Fig. |3|. 

In the present study we focus on the case of a heavy quark of mass M. Comparing 
—AE (L P ) to the stationary collisional energy loss —AE co u(L p ) oc L p [19], our esti- 
mate (see (|J) and Fig. |) (-AE coll (L p ) + AE {L P ))/E = ^/2/3C F a s m D /M ~ 5%, 
for a charm quark of mass M = 1.5GeV and for L p > 4fm, is numerically consistent 
with [13,20]. However we find explicitly that the reduction of —AEq(L p ) compared 
to —AE co ii(L p ) scales as 7 when 7 — > 00, which can be interpreted as an effective 
'retardation' of the stationary regime (see Fig. |3[) scaling similarly. We think it is 
important to extract such an asymptotic parametric behaviour in order to get a 
better insight on heavy quark energy loss, and also in view of applications to heavy 
ion collisions at higher energies such as those planned at the Large Hadron Collider. 
We recall that our calculation includes, in a unique classical framework, transition 
radiation, initial bremsstrahlung and collisional processes. The latter were not in- 
cluded in [12] and were calculated in [20] without the contribution from transition 
radiation (calculated separately in [13]). Let us mention Ref. [21], where an effect of 
interference between elastic and radiative amplitudes is studied. This effect is not 
included in our approach. 

We stress that our results should be taken at the qualitative level only. Indeed, 
we work in an ideal framework, where the medium created in the heavy ion collision 
is assumed to be a static, equilibrated, high temperature QGP with small coupling 
g <C 1. We also work in the fixed coupling approximation. In the case of a running a s 
we expect the slope (—dE/dx) co u of the stationary collisional loss to be independent 
of E when E — > 00 [22], contrary to the result for fixed a s [16,19] which is logarithmic 
in E (see ( ft.lOQ ). Despite this reduced energy dependence, {—dE j dx) co u with running 
a s might be larger than in the fixed coupling approximation, because it behaves as 
(—dE/dx) co u oc a s instead of oc ot 2 s [22,23]. We thus cannot exclude that a larger 
(—dE/dx) co u than usually assumed partly compensates the drastic suppression we 
find for —AE (L P < 5fm). A consistent calculation of — AE (L p ) (and of the slope 
of its asymptote (—dE/dx) co u) with running a s is needed to answer this question. 

2. Infinite medium: theoretical model and critical discussion 

In this section we first recall the model used in [15] (which was directly adapted from 
Ref. [16]) to study the mechanical work W undergone by a heavy quark produced 

3 For a charm quark the assumption ihad ^> L p should be reasonable in practical applications 
where L v < 5 fm and E > 10 GeV. 
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at t = in an infinite QGP, and travelling the distance L in the plasma. We then 
simplify this model in order to focus on its essential feature. 

We start from the expressions of the longitudinal and transverse (chromo-)electric 
fields in momentum space K = (u, k), 

El = — ~ ; e t = -. %n- 2 , (2-1) 

which follow from Maxwell equations in a medium with dielectric functions and 
6t, in the abelian approximation and within linear response theory. For a heavy 
quark of color charge q a (with q a q a = Cpa s ) produced at t = 0, the classical 3-vector 
current density j a is of the form [15]: 

f{K)=iq a V , , (2.2) 
K.V + if] 

with longitudinal and transverse components given by ji = (j.k/k 2 )k (where k = \k\) 
and ]t = j — Jl- The quark 4- velocity is denoted by V = (l,v) and assumed to be 
constant in the following. Let us mention that we can easily form a conserved 4- 
current whose vector component is precisely given by (|2.2j ), for instance by assuming 
the fast heavy quark to be produced in conjunction with a static antiquark at t — 0. 
This was done in Ref. [15], where it was shown that the dominant effect (scaling as 7 
when 7 — > 00, see ( |2.11|) below) on W arises from the fast quark contribution to the 
conserved current, given by ( |2.2j ). The precise form of the conserved current used 
in [15] only affects the longitudinal contribution to the medium-induced energy loss, 
which receives a term corresponding to the difference of the dipole binding energies in 
medium and in vacuum. When 7 ^> 1 this contribution is subleading [15] compared 
to the dominant (transverse) contribution oc 7. In the present study, where we focus 



on the leading effects scaling as 7, we can thus consider (|2.2|) as a relevant model for 
the current. 

In coordinate space the medium-induced electric field reads 



d A K 



S a (t, x) = I f^- e'^-^ E a L + E^ , (2.3) 



(27T) 



ind 



where the 'induced' prescription corresponds to subtracting the vacuum ei = €t = 1 
contribution. The induced mechanical work W of the electric force on the quark is 
given by 



i-l/v f d 4 K f L / v 

-W = -v- I dtq a £ a (t,vt) = -q a v- / -— - \ dte~ lKyt 







W Jo 



Ej + E% 



J ind 

(2.4) 

In Ref. [15] the work W is calculated by choosing for the dielectric functions 
and €t those of a high temperature QGP. The latter have a rich analytical structure, 
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which somewhat complicates the calculation of the energy loss. In particular, the di- 
electric functions have a cut in the spacelike \u\ < k region (physically corresponding 
to Landau damping), responsible for the leading large L behaviour, 

-W(L) ~ -AE coll (L) + O(L ) where - AE^L) = ( ~ ) L .(2.5) 

L ^°° V ^ ) coll 

We are interested in the difference, denoted as d(L), between —W(L) and the sta- 
tionary law for collisional loss —AE co u(L), 



d(L) = -W{L) + AE coU (L) = -W{L) + (*f) L 

\ax J coll 



(2.6) 



By definition the function d(L) tends to a constant when L — > oo. 

It has been shown [15] that for 7 > 1, d(L) arises dominantly from the transverse 
contribution and from the phase space region k ~ and x = u/k — > 1. In this 
limit the HTL [17] transverse gluon propagator takes a simple form, 



At(uj = kx, k) 



10 ' 



k 2 — Ht(x) 6T^J 2 — k 2 x-+i uo 2 — k 2 



(2.7) 



where m = m£,/\/2 is the asymptotic (transverse) gluon thermal mass and the 
retarded prescription to — > to + ir/ is implicit. Thus d(L) can be simply modelled by 
assuming a medium with dielectric function 



1 - 



m 



(2.8) 



All medium effects are encoded in the single parameter m. As we will shortly see, 
the model ( |2.8| ) captures the exact leading order in 7 ^> 1, and will greatly simplify 
the calculation of transition radiation in the next section, due to e(u) being indepen- 
dent of k - corresponding to the plasma response involving no spatial dispersion 4 . 
The main effects included in our calculation - initial bremsstrahlung and transition 
radiation - will arise from the domain u ~ k oc 7m, justifying a posteriori using 
Q2.8|) . We evaluate ( |2.4|) by performing the uj- integral using Cauchy's theorem, and 
the function d(L) by subtracting the leading term (linear in L) when L — > 00, 



d{L) 



(E k -k-v)L/(2v) 



(E k -k- v) 2 



sin 



(m-k- v)L/(2v) 



(m — k ■ v)' 



— S(m — k-v) 
2v V ; 



(2.9) 



ind 



The isotropy implied by (2.8) obviously imposes £l = £ t = 
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Figure 1: The function d(L) defined by Q2.6| ) and given by ( |2.9| ) in the model Q2.S| ) (thick 
lines). We consider a charm quark of mass M = 1.5 GeV, use a s = 0.2, m 2 = m 2 D /2 = 
2ira s T 2 (l + nj/6) ~ (0.32 GeV) 2 for rif = 2 flavours and a plasma temperature T = 
0.25 GeV. We show d{L) for two values of the charm quark momentum, p = 10 GeV and 
p = 20 GeV. The thin straight lines give the values of (see ( 2.10) )) for the corresponding 
values of p. The (HTL) results of Ref. [15] for d(L) are given by the dashed lines. 



where = \/k 2 + m 2 and the 'induced' prescription corresponds to subtracting 
the value of the integrand at m = 0. In ( |2.9| ) the second line corresponds to the 
longitudinal contribution, which will be subleading in the following. 
The explicit calculation of ( |2.9| ) yields: 

doc = nm d(L) = —2Cpa s m 

L^oo 

where the last term in the brackets stands for the longitudinal contribution. In the 
v — *■ 1 limit, (|2.10| ) reproduces exactly the result of [15] obtained within the HTL 
approximation, and arises dominantly from the transverse contribution 5 : 

doo ~ -V2C F a s m D j . (2.11) 

v— >1 

The numerical evaluation of d(L) is displayed in Fig. |l|, where it is also compared 
to the results of [15] obtained with the HTL expressions of the dielectric functions. 
Since the results of Ref. [15] for d(L) are very accurately reproduced with the simple 
model (|2.8|) , we will use the latter model in section ^| in order to derive the 'zeroth 
order' of the heavy quark energy loss —AEq(L p ) in a finite QGP. 

Before doing that let us discuss (see also [18]) the incorrect interpretation which 
was done in [15]. In [15] the stationary collisional loss —AE co u(L) oc L [19] was added 
to d(L), and the result interpreted as the quark energy loss —AE(L) in a infinite 

5 We note that the longitudinal contribution given in the second line of ( |2.9|) is not the same as 
in Ref. [15]. In particular it does not include the induced binding energy of the dipole created at 
t = 0. This is because instantaneous (ui = 0) self-interactions are suppressed in the model (|2.S| ). 
These details turn out to be irrelevant since the longitudinal contribution is subleading when 7 > 1, 
both in Ref. [15] and in the present study. 



.vT 



arcsm v 71 
v 2v 



(2.10) 
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plasma. Due to the large negative value of doo, a significant delay before the onset 
of the stationary (linear) behaviour was observed, and designated as 'retardation 
effect'. As noted in [15] this 'retardation effect' is not a genuine retardation of purely 
collisional loss, since various physical effects contribute to the apparent delay. First, 
part of the work done on the charge is due to the initial radiation induced by the 
sudden acceleration of the quark at t — 0. Due to the difference between the gluon 
dispersion relations in medium and in vacuum, initial bremsstrahlung in QGP differs 
from that in vacuum. This is the QCD equivalent of the Ter-Mikayelian (TM) effect, 
studied in [14] and also properly identified in [15]. In the limit 7 ^> 1, half of the 
apparent 'retardation' is actually due to the TM effect [15]. 

As can already be seen in the case of a charge with v ~ 1 produced in vacuum, 
initial radiation represents only half of the mechanical work done on the charge 
after its production. Using Poynting's theorem we can show that the other half 
is associated to the creation of the charge's proper field. Contrary to what was 
implicitly done in [15], this self-energy contribution should not be counted as 'energy 
loss' as it is part of the charge asymptotic state. With an accurate definition of 
energy loss in terms of asymptotic states [18], we find that the retardation time of 
purely collisional energy loss is quite small, t ret ~ l/m D , instead of t Tet oc •y/m D as 
argued in [15]. 

In the next section we will consider the more realistic situation, where the quark 
is produced in a finite size QGP. In this case the 'asymptotic' (we assume the quark 
to hadronize long after escaping the medium) quark self-energy is the same as in 
vacuum. Interpreting (minus) the mechanical work on the quark as energy loss is 
then legitimate, provided transition radiation is taken into account. Let us stress 
here that our final result for —AE (L P ) in Fig. § displays an effective time delay 
t re t oc 'j/m£, before the linear regime. This effective delay is not due to a retardation 
of collisional energy loss (at least at leading order in 7 ^> 1), but arises from a 
non- compensation - already noted in [13] - between the TM effect and transition 
radiation. 

3. Finite size medium: implementing transition radiation 

Here we study how the finite size L p of the medium affects the work W done on 
the heavy quark. Instead of W(L) = W(L,L P = 00), we now consider the explicit 
dependence of W(L, L p ) on L p , and implement transition radiation induced by the 
discontinuity at L = L p between medium and vacuum. 

We derive transition radiation along the lines of Ref. [24]. We assume that 
the separation surface between medium and vacuum is a plane located at z — L p , 
and that the heavy quark produced at t — and z = in the medium (with the 
associated current ( |2 .2| ) ) travels along the z-axis (with v z = v > 0). We denote the 
medium and vacuum as media 1 and 2, with dielectric functions €i(u>) given by ( |2.8j ) 
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and e 2 (uj) = 1. (This is similar to the model used in [13]). Transition radiation 
is obtained by adding to the fields (|2.1|) in both media some (transverse) fields E' T 
solutions of the homogeneous (abelian-like) Maxwell equations. It is convenient [24] 
to express these fields in the mixed (uj, k±, z) representation, 



4vrg a hi(u, k ± ) e - iax{ - z ~ L ^ (k\e z + aifcj 



E&(u), k±, z) = 4nq a h 2 (u, k ± ) e ia2iz ~ Lp) (k\e z - a 2 k^j (3.1) 



Oi = aAu, k± ) = \ €i u 2 - k 



where the subscript "_L" denotes a vector component orthogonal to the unit vector 
e z specifying the z-axis, the retarded prescription u — > uj + irj is implicit, and the 
square root is defined as having a cut along the positive real axis, i.e. such that 
Im y/z > for all complex z. One directly sees that the fields satisfy 6 the 

(homogeneous) Gauss and d'Alembert's equations, which read d z E z + ik±_ ■ E± = 
and (euJ 2 — kj_)E + d\E = in the mixed representation. 
With our definition of the square root we have 

sign(u;) uj 2 — k\ for e« uj 2 > k\ region (i) 

*i = { t • (3-2) 

i^y kj_ — eiLU 2 for ^u 2 < k\ region (ii) 

The electric field in coordinate space arises from either region (i) or (ii) in momen- 
tum space. The contributions from region (i) correspond to a superposition of plane 
waves propagating with decreasing z (for E' T ^) or increasing z (for E' T2 \ in a direction 
specified by the wave vector k = (k±, —<T\) or k = (/c±,cr 2 ). The absence of compo- 
nents cx e +laiz in E' TX and oc e~ lr72Z in E' T2 is dictated by the transition fields being 
created at the interface between medium (z < L p ) and vacuum (z > L p ), so that 
no wave comes from z = — oo in the medium or from z = +oo in the vacuum. The 
contributions from region (ii) correspond to surface waves which decay exponentially 
in the longitudinal distance. Including terms oc e +l<JlZ in E' TX and oc e~ 1<J2Z in E' T2 
would be unphysical as they diverge for z — > — oo and z — > +oo respectively. We 
stress that both types of waves (plane waves and surface waves) are simultaneously 
included in our framework by defining the square root with a cut along the positive 
real axis (with the implicit retarded prescription to — > u + irj) . 

The functions h\ and h 2 in ( ^.1| ) are obtained by requiring that the total (inho- 
mogeneous) fields 

rpa po i i?a i 7? 'a 

rpa 7?a 1 7?a 1 rp'a (o q\ 



6 Recall that the permittivity does not depend on k z in the model (2 
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are consistent with the continuity of the components E±(t, x±, z) and D z (t,x±, z) 
at the transition surface z = L p . We evaluate hi(u,k±) and h,2(w,kx) explicitly in 
Appendix A. The results are given by (|A.8| ) and ( |A.10| ). 

The work —W(L, L p ) is given by an expression similar to (|2.4| ), but we now have 
to specify whether the quark is in medium 1 or medium 2 (the vacuum) at time t: 

-W(L,L P ) = -q a v- j' dt {e(L p -vt)£?(t,m)+Q(vt- L p )£%(t,vt)} . (3.4) 

The medium-induced fields are related to the total fields ( |3~3| ) which include the 
transition fields. Instead of (12.31) we have 



£f(t,x) = I ^Le-^-^ E^ + E^ + E^ _ _ . (3.5) 



(27T) 



iud 



Inserting (|3.5|) in ( |3.4|) and subtracting the leading term (linear in L) when L — > oo 
we get: 

d(L) - W'(L, L p ) = d'(L, L p ) ; L = Min{L, L p } , (3.6) 

where d(L) is given by (|2.9| ) and the contribution —W'(L, L p ) arises from the tran- 
sition fields: 

-W'{L,L p ) = -f*. J^fj^dte^M^ 

r-L/v 



+Q(L-L p ) dte- lKVt E T %) . (3.7) 

J Lp/v 

Using the mixed representation of the transition fields and integrating over time 
we obtain: 

-W'(L,L P ) = -f-i / dk\kl / du\ C — iht(u;,k ± )e^ L * 

+e(L-L p ) i/u(w,*i)e"^ . 

(3.8) 

The function d'(L, L p ) defined in (|3.6j ) is shown in Fig. |2|. For finite L p , d'(L, L p ) 
has a cusp at L — L p , where the effect of the transition fields on the work done on 
the heavy quark starts 7 . Most importantly, the asymptote d'^ of d'(L, L p ) in the 
limit L ^> L p — > oo scales in 7 when 7 — > 00 (see (B.8)), 

d 'oo - -doo ^ --V2C F a s m D >y . (3.9) 
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Figure 2: The function d'(L,L p ) of (gl| for L p = 2, 4, 6, 8 and 10 fm (from top to 
bottom) together with the function d(L) (lower curve) of Fig. [j], for p = 20 GeV. The 
dashed straight lines represent the asymptote d'^ of d'(L, L p ) for the different values of L p . 



As discussed in the end of section in the case of a finite plasma size it is 
legitimate to interpret the mechanical work d'{L = oo, L p ) as quark energy loss. 
It is also clear that (f (oo, L p ) includes the (radiative) contributions from initial 
bremsstrahlung and transition radiation. In order to obtain the 'zeroth order' loss 
—AE(L = oo, L p ) = —AE (L P ), we must add to d'(oo,L p ) the stationary law 
—AE co u(Lp) for purely collisional loss. For the latter we take the HTL result of [19] 
in the logarithmic accuracy log (/c max /m 9 ) 3> 1, 



> 

2 



1 1-w 2 , fl+V 

o - ! - lo S 1 

v 2f \ 1 — v 



log ( ^ ) , (3.10) 



m g 



— — , V ET > and m g — mo/ v^3 



The resulting energy loss 

—AE(L = oo, L p ) = -AE coll (L p ) + d'(L = oo, L p ) 

= -AE cM {Lp) + d{Lp) - W\L = oo, L p ) (3.11) 

is determined by (|2.9|) , (|3.10|) and (|3.8|) , with the functions hi and h 2 given by (|A.8| ) 



and ( WM ). 

Our result for —AE (L p ) is shown in Fig ^. We observe an important sup- 
pression of the zeroth order energy loss as compared to the stationary collisional loss 
— AE co ii(L p ) . This is mainly due to the fact that the asymptote d'^ of d'[L = oo, L p ), 
though reduced by a factor 1/3 compared to (see ( |3.9| )), still scales as 7 when 
7 ^> 1. As already mentioned in the end of section 0, this is due to a non- 
compensation between the TM effect and transition radiation. As is obvious from 



7 Strictly speaking, the effect of transition on d'(L,L p ) starts at L = L s = 2vL p /(l + v) < L p , 



in agreement with causality, as shown in Appendix B, sec (B.l). For p > 10 GeV, v ~ 1 and the 
difference between L s and L p cannot be seen on Fig. ^|. 
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Figure 3: Charm quark energy loss —AEq(L p ) = —AE(L = oo, L p ) for p = 20 GeV (solid 
line) given by ( 3.11| ). —AEq(L p ) is compared to the stationary law ( |3.10| ) for collisional 
energy loss —AE co u(L p ) (straight line) [19]. 



Fig |3|, we can define an effective retardation time t mt before the linear regime for 
—AE (L p ) sets in, 

tiek = df 0O /(dE/dx) coU , (3.12) 

leading to the scaling t Tet oc r y/m D , thus to a large effective delay. Numerically, 
t Tet ~ 5 fm for p = 20 GeV. Fig [3] is our main result, which shows that invoking only 
the heavy quark collisional energy loss —AE co u(L p ) to explain the 'single electron 
puzzle' is not satisfactory. Indeed, adding the other relevant contributions (initial 
bremsstrahlung and transition radiation) to obtain the energy loss —AE (L P ) to 
zeroth order in an opacity expansion [6], we find that —AE (L P ) <C —AE co u(L p ) for 
in-medium quark path length L p < 5 — 7 fm. 



4. Conclusion 

In this paper we studied the energy loss — AE (L p ) of a fast heavy quark produced 
in a finite size QGP, to zeroth order in an opacity expansion. —AEq(L p ) includes 
collisional energy loss, initial bremsstrahlung arising from the quark creation at t — 0, 
and transition radiation appearing when the quark passes the discontinuity between 
medium and vacuum. Only the radiative loss induced by rescatterings in the plasma 
—AE ra d(L p ) [5,6] should be added to —AE (L P ) to obtain the total heavy quark 
energy loss. Our main result is that —AE (L P ) is strongly reduced compared to the 
stationary linear law for collisional loss, as seen on Fig [| Due to the negative shift 
d'ac scaling as 7, the linear regime for —AE (L p ) sets in after a quite large effective 
retardation time t ret oc 7. 

From Fig. [| we conclude that for the generic values L p ~ 5 fm and p = 20 GeV, 
the loss —AE (L P ) should be negligible, and the total heavy quark energy loss should 
be theoretically correctly estimated by —AE rad (L p ) only. Our results suggest that a 
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missing contribution to the heavy quark energy loss, as proposed in Ref. [11], might 
not be a correct explanation for the surprisingly strong nuclear attenuation of large 
p± electron spectra from heavy flavour decays. 

We already mentioned in the Introduction that our results should be considered 
at the qualitative level, in particular because the slope (—dE/dx) co ii of the stationary 
regime usually considered [16, 19] is modified when the running of a s is taken into 
account [22]. As discussed in [15], the small L (small L p ) behaviour of Fig. [T] (Fig. £J) 
might also be affected by the running of a s . However the position and scaling in 
7 of the asymptotes d^ and d'^ (see ( |3.9|) and Fig. |3|) should be unaffected by the 
running of the coupling [15], leaving unchanged the main result of the present study. 
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A. Derivation of the transition fields 

In this Appendix we obtain the functions h\(u, k±) and h 2 (u, k±) defining the tran- 
sition fields (|3TT|). Written in the mixed space (a>, k±, z), the continuity conditions at 



E 2 n (co, k±, L p ) + Et 2 (u, k±, L P )J (A.l 
D 2 n (u, k±, L p ) + e 2 E^ 2 (u, k ± , L p ) (A.2 



the surface z = L p read: 



E™(u, k±, L p ) + Etx(uj, k ± , L p ) 
D™(uj, k±, L p ) + exE' Tl (to, k±, L p 



where the color index is implicit and the inhomogeneous fields are given by (i = 1, 2) 

dk z 



E™(u,k ± ,z) = I —e 



E u (u, k±, k z ) + E Ti (uj, k±, k z 



D] n (u,k ± ,z) = ei(u)El n (uj, k±, z) 



(A.3) 
(A.4) 



Using the definition (|3.1|) , we can easily solve for the functions hi and h 2 . We find 
from the above continuity conditions: 

£2 - ei 



with 



J 



hi(u>, k±) 
h 2 (uj, k±) 

q a k± ■ \E 2 na (L p ) - Ef a (L % 



e\fJ2 + 620"! 

fi - ^2 
ei(T2 + 62(71 



[e 2 J + a 2 K] 
[-e 1 J + a 1 K] 



(A.5) 



ArKC F OL s k\{e 2 - ei) 



and K 



Df a {L p ) - D [ r(L p 



47rC F a s k 2 , (e 2 - e x ) 



(A.6) 
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Using fl2.1|) (with eu = en = q) and (|2.2|) we obtain: 

dk z e ikZL " k z [(ei + e 2 )u 2 - A; 2 ] 



J(uj, k\ 



v 



ue\e2 J Itx (uj — vk z )(eioo 2 — k 2 )(e 2 uj 2 — k 2 ) 
dk z e *k z L p 

~'~ ' 2vr (u - vk z )(enu 2 - k 2 )[e 2 uj 2 - k 2 ) 



(A.7) 



From ([A. 5 ) and ( |A.7| ) we readily check that hi(—uj,k±) = hi(u,h±)*. Hence the 
transition fields ( |3.1| ) satisfy E^(—uj, —k±, z) = E^(cu, k±, z)*, ensuring the reality of 
the fields in coordinate space. We can further put (|A.5|) in the form 



hl\LL),k±) = l ; li 



ho(u, k 



e\o 2 + t 2 o\ 
(ei - e 2 )uj 



i- 



e 1 a 2 + e 2 0"i 
where the functions Ii and I 2 read 

dk z e tk z L p 



(A. 



h = 

1 2 = 



eiur 

iv 
e 2 uj 2 



2n (to — vk z )(e\UJ 2 — k 2 ) 
dk z e ikZL " 
2tt (u - vk z )(e 2 uj 2 - k 2 ) _ 



k z {k z + a 2 )-e 1 iu 2 

k z + a 2 
k z (oi - k z ) + e 2 uj 2 

k z — o\ 



(A.9) 



We calculate l\ and l 2 by closing the contour of the /^-integral in the upper half of 
the complex plane. We obtain: 



h 



e^(^-e 1 uj 2 + ^) 



1 

e 2 u; 5 



+ 



e i ^ L vt ( y 1( y 2 - k 2 . 



[£ - °l] K + **] 2(71 h " S fa + <*0 

^ L p( ai a 2 + h\ 



e l ^(a^ + e 2 LU 2 -^) 



+ 



£2e 



at 



] [s - ax] 2a 2 [tr 2 - a] ( ff2 - ( Cl - e 2 ) [a x - s] 

(A.10) 



The transition fields (|3.1|) are determined by ( |A.8| ) and (|A.1Q ), where we recall that 
the retarded prescription u — > a; + zt? is implicit, and that the square root satisfies 
Im sfz > for complex z. 



B. Some properties of —W'(L,L p ) 

Here we study some features of the contribution —W'(L, L p ) to the heavy quark 
energy loss (see (|3.8|) and ( |3.11| )) which arises from the transition fields ( |3.1| ). We 
first show that this contribution is consistent with causality. Then we derive, for 
7 ^> 1, the limit of —W'(L,L P ) for L ^> L p — > oo. This limit corresponds to the 
amount to be added to the asymptotic value d^ of d(L) to reach the asymptote d'^ 
of d'(L) (see Fig. §. 
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Causality 



From causality the heavy quark cannot be sensitive to the transition fields right 
after its production at t — 0. The effect of transition fields can start only when a 
signal emitted by the quark at t = (and z = 0) along the positive z-axis has had 
time to reach the medium- vacuum separation surface located at z = L p , and to come 
back to the position of the quark. Since the quark is moving with velocity v, and 
the largest speed of light in a medium of permittivity ( |2.8| ) is unity (corresponding 
to uj — > oo), this time is given by t s = L s /v = 2L P — L s . Thus —W'(L,L P ) must 
vanish when the quark has travelled less than the distance L s , 

L<L s = -^^-W'(L,L p ) = . (B.l) 
1 + v 



The latter causality requirement can be proved directly from ( |3.8|) . For L < L p 
only the first term in the bracket of ( |3.8| ) contributes. Inserting the expression of hi 
given by ( |A.8|) and ( |A.10| ) we obtain: 



Ct?o f°° f°° du p i( ~ii+ a ^ L p - e i (^+^)( L p- L ) 

-W>(L,L P )\ L<L =^ dklkl 1 <U< 



2tt Jo J-00 ^ f + °"i 



eia 2 + e 2 a t \ [f - <j\] [a + a 2 ] 2a x [a t - a] {<r x + tr 2 ) < 

Recalling (see section 3) that the square root satisfies Im \fz > for complex z, we 
can perform the k>integral by closing the integration contour in the upper half-plane 
provided the phase 

2ai L p - (- + <n) L (B.3) 

has a positive imaginary part when |a>| — > oo in this half-plane. Since for Imw > 
the integrand of ( |B.2| ) has no singularity (which can be easily checked) and <j\ ~ uj 
when |o;| — > oo, we obtain that —W'(L,L P ) vanishes when 2L P > (1 + l/v)L, which 
proves ( |B.1| ). 

Limit of —W'(L, L p ) when L L p — > oo 

We first consider the limit of —W'(L, L p ) for L — > oo at fixed L p . This represents 
the amount brought by transition radiation, to be added to — AE coll (L p ) + d(L p ) (see 
( |3.11| )) in order to get the heavy quark energy loss —AE (L p ) represented in Fig. [| 



For a fast quark, 7 > 1, this limit is given by the second term of ( |3.8| ), where the 
term with a rapidly oscillating phase factor oc e~ l( -« _cr2 ^ L is neglected: 

-W\L 00, L p ) = %i / dkl kl [ du Tu e V _*? 2 ~ e , l)U _, h , (B.4) 



2tt 7 J~oo " [f - 0-2] (ei<7 2 + e 2 «7i) 
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where we used ( |A.8 ) and I2 is given in (|A.1Q ) 



We now evaluate the limit of ( |B.4| ) when L p — > 00. In this limit only the first 
term of I2 in ( |A.10| ) contributes to (|B.4|) , the two other terms oscillating rapidly. We 

get: 



-W'(L,L P ) 



2tt 



L»L p -+oo 

V 



i-00 W ( e , ff2 + e . 2(Tl ) [K + j [H - (7 2 J [J - ff J 



The latter integral can be evaluated for 7 ^> 1 by anticipating that it arises domi- 
nantly from the domain u ~ 7m and k± ~ m. Approximating the integrand in this 
region and using (|2.8| ) (and €2 = 1) we obtain 



L>L,-aa n J Q ± ± J ^ + fc 2 )2( £ + fc 2 + ^ 

(B.6) 

The latter double integral can be exactly evaluated, yielding (use m = m D /\/2): 

-W'(L, L p ) \ L>>Lp ^ = 4, - ^ ^ I y/7C F a.m J r t • (B.7) 

The expression ( |B.6| ) is also trivially shown to be dominated by u ~ 7m, fc^ ~ m, 
justifying our initial approximation. The result ( |B.7| ) is exactly 2/3 of \doo\, see 
( CT) . Hence: 

C - o^oo ^ --V2C F a s m D "/ . (B.8) 



Finally, we note that in ( |A.10| ) the second and third terms of I2, which have been 



neglected to extract the large L p limit of (p3.4[) , are proportional to the phase factors 



e H„ <ri,2)L p _ w^^gn ^ j m ^> m anc i ~ m we can approximate, for instance, 



,2 



- (Tij L p ~^\^— + ki + m 2 j~0 {L v mli) , (B.9) 

Thus the large L p limit used above should be understood as L p ^> j/tud, where the 
neglected terms are indeed rapidly oscillating. In this limit the function d'(L, L p ) 
is close to its asymptote d'^, i.e., the energy loss —AE (L p ) including the effect of 
transition radiation (see Fig. |3|) is close to its asymptotic regime. 
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